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Abstract 

We propose an algorithmic procedure i) to study the "distance" 
between an integrable PDE and any discretization of it, in the smaU 
lattice spacing e regime, and, at the same time, ii) to test the (asymp- 
totic) integrability properties of such discretization. This method 
should provide, in particular, useful and concrete informations on 
how good is any numerical scheme used to integrate a given integrable 
PDE. The procedure, illustrated on a fairly general 10-parameter fam- 
ily of discretizations of the nonlinear Schrodinger equation, consists of 
the following three steps: i) the construction of the continuous mul- 
tiscale expansion of a generic solution of the discrete system at all 
orders in e, following ii) the application, to such expansion, of the 
Degasperis - Procesi (DP) integrability test [21 [3], to test the asymp- 
totic integrability properties of the discrete system and its "distance" 
from its continuous limit; iii) the use of the main output of the DP test 
to construct infinitely many approximate symmetries and constants of 
motion of the discrete system, through novel and simple formulas. 

1 Introduction 

Given a partial differential equation (PDE) and a partial difference equation 
(PAE) discretizing it, it is interesting to know, when the lattice spacing 
e is small, "how close" the two models are. In particular, if the PDE is 



1 



integrable, it is important to have a way to establish if such a discretization 
preserves integrabihty or, at least, how "close" is to an integrable system, 
detecting the order, in e, at which the discretization departs from integrabihty 
and, correspondingly, the time scale at which one should expect numerical 
evidence of nonintegrability and/or chaos. In addition, given a PDE and two 
PAEs discretizing it, it is also interesting to know, when the lattice spacing 
e is small, "how close" the two PAEs are. 

In this paper we propose to answear these basic questions in the following 
way. Concentrating on an integrable PDE and on a PAE discretizing it, 

1) we construct and study in detail the multiscale expansion at all orders 
of a generic solution of the PAE under scrutiny, generated in the small 
e regime, following the procedure developed in [I]. At 0(1), the leading 
term u of such asymptotic expansion satisfies the integrable PDE; to keep 
the expansion asymptotic, we eliminate the secularities due to the linear 
part of the PAE, arising at each order, introducing infinitely many slow 
(time) variables and establishing that the evolution of u with respect to such 
slow times is described by the infinite hierarchy of commuting flows of the 
integrable PDE, as in [T]. 

2) We make use of the asymptotic integrabihty test developed by Degasperis- 
Procesi (DP) in [21 [3] on such a multiscale expansion to test, at all orders, the 
"asymptotic" integrabihty properties of the PAE; in particular, detecting the 
order in e (and, correspondingly, the time scale) at which the discretization 
departs from integrabihty. At this time scale, f.i., numerical simulations are 
expected to give some evidence of non- integrable and/or chaotic behaviour. 

3) We flnally show how to make use of the main output of the DP test to 
construct inflnitely many "approximate" symmetries, at a required order in 
e, of the PAE under scrutiny, using novel and simple formulas. 

Recent studies on the performances, as numerical schemes for their contin- 
uous limits, of PAEs possessing the same (continuous) Lie point symmetries 
as their continuous limits can be found in [11[5]. Studies on the performances, 
as numerical schemes for their continuous limits, of integrable discretizations 
of integrable PDEs can be found, f. i., in [B] and [7|; in this case, the inte- 
grable discretization possesses inflnitely many exact generalized symmetries 
and constants of motion in involution at any order in e, reducing to the gen- 
eralized symmetries and constants of motion of the integrable PDE in the 
continuous limit. The PAEs selected by our approach possess instead in- 
flnitely many approximate generalized symmetries and constants of motion 
in involution at the required order in e (see §3.1), reducing to the generalized 
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symmetries and constants of motion of the integrable PDE in the continuous 
hmit. 

The procedure we propose should allow one to have a control on the "dis- 
tance" between the PAE and its continuous limit, as well as on the distance 
between two different discretizations of the same PDE. Indeed, suppose we 
construct an asymptotic expansion of the form ip = u + 0(e"), a > 0, where 
is a. generic solution of the PAE and u is the corresponding solution of its 
continuous limit; if, at O(e^), /? > 0, the PAE passes the DP test, we infer 
that II'?/' — m|| = 0{e°') at time scales of 0(e~^), where || ■ || is the uniform 
norm wrt x and t (the norm used to test the asymptotic character of the 
generated multiscale expansion). In this way, since we control the distance 
between "generic solutions" of the PAE and of its continuous limit, we also 
control the distance between the PAE and its continuous limit. In addi- 
tion, if the multiscale expansions of two different discretizations of the same 
PDE pass the DP test at O(e^), we infer, from the triangular inequality, that 
1 1?/' — 01 1 < 1 1-?/" — m| I + 1 10 — n| I = 0(e°) at time scales of 0(e~^), where ip, cj) 
are solutions of the two different discretizations of the PDE corresponding 
to the same generic initial-boundary data; therefore we have a control also 
on the distance between the two different discretizations of the same PDE. 

Some historical remarks are important, at this point, on the theory of 
multiscale expansions in connection with integrable systems, to put the re- 
sults of this paper into a proper perspective. Multiscale expansions of a given 
PDE are very useful tools for investigating the properties of such a PDE and 
for identifying important model (universal) equations of physical phenom- 
ena. For instance, if the original nonlinear PDE has a dispersive linear part, 
a small amplitude monochromatic wave evolving according to it develops a 
slow space-time amplitude modulation described by the celebrated nonlinear 
Schr5dinger (NLS) equation [HI El [IDl [IB [I2| (see also [SI El [IS]) 

-|- -|- 2c|upM = 0, u = M(x,t)GC, (1) 

integrable if c is a real constant [16j. Considering, instead, three monochro- 
matic waves and imposing a suitable resonance condition on their wave 
numbers and dispersion relations, one generates another integrable univer- 
sal model, the 3- wave resonant system [17J. In the above two examples, 
the expansion is constructed around "approximate" particular solutions of 
the original PDE (the monochromatic waves). It is also possible to expand 
around "exact" particular solutions of the original PDE; for instance, as 
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shown in [IB], expanding around the exact solution uq = exp{2ict) of ([T]), 
the first nontrivial term of the asymptotic expansion evolves according to 
another important model equation: the Korteweg-de Vries (KdV) equation 
[19j, sharing with NLS the property of integrability [20j. Since multiscale ex- 
pansions preserve integrabilty [18], i) if the original PDE is a "C-integrable" 
system (i.e., it is linearized by a "change of variables" [13121], like the Burg- 
ers equations [22]), the model equation generated by it is linear [T51 [21]; ii) 
if the original PDE is an "S-integrable" system, or soliton equation (like the 
NLS equation), integrated in a more complicated way via a Riemann-Hilbert 
or ^-problem [231 1211 [2S1 126] , the model equation generated by it is also "S- 
integrable" ; viceversa, iii) if the model equation generated by the expansion is 
not integrable, then the original equation is not integrable too. This criterion 
has been used in [151 123 1211 129] as a simple test of integrability. In addi- 
tion, the universal character of the identified model equations (NLS, KdV 
or others) is also the reason why model equations possess very distinguished 
mathematical properties and, often, they are integrable [151 123 121] ■ 

Multiscale expansions can also be carried, in principle, to all orders and, 
as a consequence of eliminating the secular terms at each order, a sequence 
of slow time variables tn = must be introduced and the dependence of 
the leading term of the expansion on such slow times is described by the 
hierarchy of commuting flows of the integrable model equation [1]. This 
multiscale expansion at all orders has been used in [21 E] to build an efficient 
asymptotic integrability test for the original PDE (see §3 for more details 
on such test). An alternative asymptotic integrability test, based on the 
existence of approximate symmetries for the original PDE, can be found in 
[30] . The ideas and procedures developed in [H [21 [3] have been recently used 
to build an integrability test also for PAEs [311 [321 [S3]; in this approach, 
one expands, as for the PDE case, around approximate or exact particular 
solutions of the PAE under investigation, obtaining a continuous multiscale 
expansion at all orders, following [1], and applying on it the DP test. The 
main difference between the procedure followed in [311 [321 E2] and the results 
of this paper is the following. The standard multiscale approach used in [211 
[321 133] , obtained expanding around approximate or exact particular solutions 
of the PAE under investigation, cannot give informations on how close this 
PAE and its continuous limit are, the main goal of the present paper. The 
common features of the procedure in [311 1321 133] and of that used in this 
paper are that, in both cases, one constructs, from the given PAE, continuous 
multiscale expansions carried to all orders, as in [T] , and one applies to them 
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the DP integrability test. Therefore both procedures can be used to test the 
integrabihty and the asymptotic integrabihty of the original PAE. 

Another integrabihty test for PAEs is the so-caUed "symmetry approach" 
[34j , based on the existence of higher order symmetries and originally devel- 
oped to test the integrability of PDEs [351 136] . 

The results of this paper are illustrated on the basic prototype example 
of the NLS equation ([T]), starting from the following discretization of it: 

#n,t + e"^ + " 2V^„) + F(V'n-l, ^n, ^n+l) = 0, 

F{lljn-l,i^n,'4'n+l) ■= 2ai\lpri\^1pn + a2\lpn? {lpn+l + 1pn-l) + 
as^jlilpn+l + -ipn-l) + a4j/'n(|^n+lP_+ \lpn-l\'^)+ (2) 
a5^n(^n+lV'n-l + ^/'n+lV'n-l) + ^eV'nlV'n+l + V'n-l) + 
2a7?/'„Vn+lV'n-l + asG^n+lPV'n+l + | V'n-1 PV'n-l) + 
aglV'n+lV'n-l + ^l_illJn+l) + CHO ( I + |^n-lP^n+l), 

where the constant coefficients aj, j = 1, ... ,10 are real, reducing to (JT]) in 
the natural continuous limit in which the lattice spacing e — >^ and ne — > 

a; G M, ipn{t) — > u{x,t), with 

10 

c = ^aj. (3) 

The 10-parameter family of equations ([2]) has been recently taken in [37] 
as the starting point of an analysis devoted to the identification of discretiza- 
tions of NLS that possess, at the same time, a solitary wave and a breather 
solution reducing, respectively, to the one soliton and breather solutions of 
the NLS equation ([1]), in the continuous limit e — *■ 0. We remark that, rescal- 
ing the dependent variable, one can always introduce one normalization for 
the 10 coefficients; for instance, one can choose one of these coefficients, say 
aj, to be sign{aj) or, better for our purposes, one can normalize the sum ([3]) 
of the 10 coefficients to coincide with the prescribed coefficient c of the NLS 
equation ([1]). 

The linear part of the discrete NLS (dNLS) ([2]) is the standard discretiza- 
tion of {iut + Uxx)', its nonlinear part is uniquely fixed by the following, 
physically sound, properties [37]. a) Equation ([2]) must possess the gauge 
symmetry of first kind (i.e., iiipn is a solution, il>ne~^^ is a solution too, where 
6 is an arbitrary real parameter), corresponding to the infinitesimal gauge 
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symmetry —iipn- b) The nonlinearity is cubic; i.e., is the weakest nonhn- 
earity compatible with the above gauge symmetry, c) Only first neighbours 
interactions are considered, d) Equation ([2]) is invariant under the symmetry 
transformation ipn±i if^n^i (space isotropy). 
The dNLS ([2]) contains, in particular: 

1) the integrable Ablowitz - Ladik (AL) equation [3H] 

i^n,t + ii'n+l + i'n-1 " 2^/'„) + a2|^„P(V'n+l + ^n~l) = 0, (4) 

for ttj = a2Sj2, j = 1, . . . , 10; 

2) the discretization 

+ e"^ ii'n+l + i'n-l " 2V^„) + 2ai\lpn\'^1pn = 0, (5) 

for ttj = ai6ji, j = 1, . . . , 10, relevant in several applications [39l HQl SH SSI 
Us] , whose nonintegrability has been recently shown in [311 [32] using the DP 
test; 

3) the discretization corresponding to 

Oio = 0-8, ai = 04 = as = oe = a? = ag = 0, (6) 

with a2,a3,as arbitrary, possessing a solitary wave as well as a breather 
solution reducing, respectively, to the one soliton and breather solutions of 
the NLS equation in the limit e ^ [37j; 

4) the discretization corresponding to 

^8 = ^3; 0,2 = 2a3, a4 = 2a6, = = ag = aio = 0, (7) 

where 01,03,00 are given in terms of physical quantities, describing coupled 
optical waveguides embedded in a material with Kerr nonlinearities [S]; 

5) the discretization corresponding to 

04 = 02, oi = 03 = Oe = 08 = 02/2, 05 = 07 = 09 = oio = (8) 

(a particular case of ([T])), appearing in the modelling of the Fermi-Pasta-Ulam 
problem |45j . 

For special values of the coefficients Oj's the dNLS equation ([2|) is Hamil- 
tonian. For instance, equations ([1]), ©, ([7]) and ([8|) are Hamiltonian [ilj . 

If < e << 1, the discrete scheme ([2|) approximates the NLS equation 
([ID,© with an error of O(e^). To study more precisely how close equations 
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([2D and ([T]) are and, in particular, the integrability properties of ([2]), in this 
paper we follow the procedure indicated in the first part of this introduction, 
obtaining the following results. 

1) Due to the structure of the vector field in ([2]), the generated e-expansion 
contains only even powers. At O(e^), the dNLS ([2]) passes the DP test iff the 
10 coefficients satisfy the elegant quadratic constraint 

10 

(ai — 3a3 — 2a4 — 605 — Sae + 807 — Sag — ISog — aio) ctj) = 0, (9) 

10 

factorized into two linear constraints. If the first constraint ~ ^ 

satisfied, we are in the C-integrability framework and the dNLS ([2]) approxi- 
mates the linear Schrodinger equation with an error of O(e^), for time scales 
of 0(e~^). If, instead, the second constraint is satisfied: 

ai — 803 — 2a4 — 605 — Sag + 807 — Sag — ISog — Oio = 0, (10) 

we are in the S-integrability framework and the dNLS ([2]) approximates the 
NLS equation ([I]),© with an error of O(e^), for time scales of 0(e~^). 

We remark that, among the ten single dNLS equations obtained choosing 
only one of the ten coefficients different from zero in ([2]), only the AL equation 
@j satisfies the constraint ([9]) and passes the test at O(e^). 

2) At 0(e'^) we have the following two scenarios. In the C-integrability 
framework, the dNLS ([2]) approximates, with an error of O(e^), the linear 
Schrodinger equation for time scales of 0(e~^) iff the following four linear 
constraints 

10 

O'j =0, ai + a2 + aQ + = 0, 04 — 05 + 2a8 — 209 = 0, , , 
j=i (llj 
a2 + 2(03 + 805 + 8a6 — 07 + as + ^ag) = 

are satisfied by the coefficients. Since one of the real a^-'s can always be fixed 
rescaling the dependent variable ip, equations ffTTj) characterize a 5-parameter 
family of discrete NLS equations ([2]) passing the test at such a high order. 

In the S-integrability framework, the dNLS ([2]) approximates, with an 
error of O(e^), the NLS equation ([I]),© for time scales of 0(e~^), iff the 
coefficients satisfy, together with the linear constraint flTU]) . the five quadratic 
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constraints (135]) . (IMI) - (1351) . Since these five constraints do not contain the 
term (02)^, they are trivially satisfied by the integrable AL equation (jl]), as 
it has to be. In general we do not expect a parametrization of such constraints 
in terms of elementary functions; however we have been able to construct the 
following two explicit examples of dNLS equations 

ii>n,t + e"^ {lpn+l + Ipn-l " 2?/'„) + fle ( " Sl'lpnl'^'ipn + ||^nP(^n+l + ^n-l) + 



(12) 
(13) 



#n,t + e"^ (V'n+l + " 2V'„) + Og - 48|?/'„PV„ " SV'n ( | V'n+i P + 

l^^-lH - 8V^„(V^„+iV^„_i + iJn+li^n-l) + lOtpnii^l+l + + 

+ i^l-li'n+l) + 6( | | ^^^^-l + |V^„-1 1 ^n+l)) = 0, 

satisfying such complicated quadratic constraints, corresponding to particu- 
lar cases in which the associated five quadrics degenerate into hyperplanes. 

These two distinguished models, passing the test at such a high order 
through the above degeneration mechanism, are obviously good candidates to 
be S-integrable discretizations of NLS. A detailed study of their performances 
as numerical schemes for NLS, and of their possible integrability structure 
(Lax pair, etc ..) is postponed to a subsequent paper. 

To obtain the above results, it is essential to use well-known integrability 
properties of equation ([T]) ( shared by all integrable systems; see, f.i. P5]-|15]) 
that we summarize here, for completeness. 

The NLS equation belongs to a hierarchy of infinitely many commuting 
fiows: 

ut^ = KM. neN (14) 

i.e., such that 

[KM,KM]l := K'MIKM] - K'^{u)[KM] = 0, n,meN, (15) 
where 

KrXu) [f] = lim — (16) 

is the usual Frechet derivative of Kn{u) wrt u in the direction /. The com- 
muting vector fields {Kn}n&n are arbitrary linear combinations, with con- 
stant coefficients, of the following basic symmetries generated by 
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the recursion relation 

Rf ■.= i{U + 2cud-\uf + uf)), ^ ^ 

where R is the recursion operator of the NLS hierarchy [50]. The basic 
symmetries used in this paper are: 

(To = —iu, (72 = i{uxx + 2c|-up'u), 

(Tg i(^xxxxxx ~l~ '2c{u Uxxxx ~l~ 6|tt| Uxxxx ~l~ ^UUxUxxx ~\~ ^UUxUxxx~\~ (-^^) 

15uUa;Ua;a;^ + + lOw^M^^) + 1 Oc^ ( 2^^ | u j^M^^ + 2mm^^+ 

5|UxpMa;a; + ^u\'^Uxx + ^i^^i^ + 6u|up|u2:p + 7u|upM^) + 20c^M|M|^j , 

and the NLS equation ([T]) corresponds to the flow Ut^ = K2{u) = (r2{u). 

Equivalently, the basic symmetries {cnjnGN are elements of the kernel of 
the "linearized" nth flow operator M„, G N, defined by 

Mnf ■■= - Kiu)[f]; (19) 

i.e.: 

M^c^m = 0, n,meN. (20) 
Due to (ITB]1 . these linearized operators commute: 

M„M^ = MmMn, n, m G N. (21) 
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The linearized operators used in this paper are: 



M2/ := ft, + 2c(mV + 2|n|V)), 

Mif := fu - Y2 f^^^^ + 2c(m(6u^./^. + 4m/^^) + u{2uj^ + 2u^f^+ 
ufxx) + (6c|m| V + + Auu^^)f + (9c|m|^ + 2|mx.P + 4mm^^+ 

4:UUxx fxx ~l~ Ux{^5Ux fxx ~l~ ^Uxxfx ~l~ '^'^xxfx '^ufxxx) ~l~ (SMx^xx"!" (^^) 
3MMxa;a;)/x)) ~l" 11(70011 Uxfx ~\~ Hl^xx/zx ~l" H/xxl^xa; ~l" ^1^xfxxx~^ 
^Uxfxxx +_QUxxxfx + Qufxxxx) + M^(5cm(6/^M^ + G/a^Ma: + 5ufxx) + 
fxxxx) + /(SOc^llil V + 10cii2(3|lXa.P + 5liMa;x) + 10cit^1ia;x + 

5(21/2^ + Sw^M^^a;) + u{70cuul + Guxxxx)) + /(40c>|6 + 35cu^ul+ 

ll\Uxx\'^ + I5cu'^{ul + 2uUxx) + QUxUxxx + AUxUxxx + QuUxxxx + 
2u{^Cu{Q\Ux\'^ + ^UUxx) + Mxxxx))) • 

At last, if c = 0, equations ([1]) and f|T71) lead to the linear Schrodinger (LS) 
equation 

IMi + Uxx = (23) 

and to its (trivial) symmetries (— i"^^9"ix). 

The paper is organized as follows. In §2 we construct the multiscale ex- 
pansion, in the small e regime, of a generic solution of ([2]), establishing, in 
particular, that the leading term of such expansion evolves wrt the infinitely 
many "even" time variables t2k '■= e^^''~^h, k G N+ according to the even 
flows of the NLS hierarchy. In §3, after summarizing the DP test and after 
showing how to use the main output of this test to construct infinitely many 
approximate symmetries of the original PAE through novel and simple for- 
mulas, we apply the DP test to the PAE ([2]), isolating the constraints on the 
coefficients a/s, j = 1, .., 10 allowing one to pass the test at time scales of 
0(e~2) and of O^e'"^), in both scenarios of C- and S- integrability. In §4 we 
summarize the results of the paper and we discuss the research perspectives 
opened by this work. In the Appendix (§5) we display the long outputs of the 
DP test, obtained using the algebraic manipulation program of Mathematica. 
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2 Multiscale expansion in the small lattice 
spacing regime 

If the lattice spacing e is small: < e << 1, as consequence of the invariance 
of ([2]) under the transformation ipn±i 'ipn^i and of the well-known formula 

= E ^^'^'/' (24) 

k=0 

only even x-derivatives appear at all (even) orders in e, implying that also the 
asymptotic expansion of ipn contains only even powers of e. Consequently, to 
eliminate the secularities appearing at all even orders in e, the coefficients of 
such expansion must depend on infinitely many "even" slow times fl]: 

r= (t2,t4,t6,---), hk:=e^^'"^h, keN+, (25) 

implying that 

dt dt, + e^dt, + e^dt, + .... (26) 

Therefore we are lead to the following ansatz for the asymptotic expansion 
of the "generic" solution of ([2]) 

^„(t) = ^e''=M('^+^)(x,t), M«(x,t) = M(x,t). (27) 

fc=0 

Plugging ( l24l) . (1261) and ( |2711 into equation ([2]) and equating to zero the 
coefficients of all powers in e, we obtain the following results. 

At the leading 0(1), we obtain the NLS equation for the leading term 
u*^^-* = u wrt the first time t2 = t: 

Ut2 = K2{u), 

10 /2g) 
K2{u) := a2{u) = i{u^^ + 2c\u\'^u), c=Y,(^j- 

j=i 

As usual in perturbation theory, at the next relevant order (O(e^) in our 
case), the "linearization" M2U^^^ of {ut^ — K2{u)) appears, together with the 
linear term {ut^ — {i / 12)uxxxx) ■, coming from the linear part of ([2]), and with 
a nonlinear term Gs, coming from the nonlinear part of ([2]): 

M2M^^^ = - [uu - i^Uxxxx^ + Gs, (29) 

11 



where 



51 = as + a4 + + as + ag + aio, ,„„v 

52 = (22 + ^4 + ^5 + 2(a6 + a? + ^8 + ctg + aio), 

53 = 2(04 — as + 2as — 2ag), S4 = 2{aQ — aj + as + ag — oio). 

Concentrating on the hnear terms in round bracket, we observe that Ut^ G 
KerM2 and {—{i/12)uxxxx) is the hnear part of the symmetry (2/A\)a4{u) G 
KerM2- Therefore, adding and subtracting the symmetry (2/4!)(T4, equation 
( !29|) is conveniently rearranged in the following way, isolating the resonant 
terms in round bracket: 



M2w(') = -K + |^4H)+^75, (31) 

where 



5f5 := i{ci\u\'^u + C2uu^^ + Csu\uJ\^ + C4\u\'^Uxx + c^v^u^x) (32) 

and 

Cl = -\c^, 

C2 = —\{ai + a2 + as + 04 + as — Sae + Say — Sag — Sag + 5aio), 

C3 = — ■|(ai + a2 + as — 5a4 + 7as + ae + ay — llag + ISag + aio), (33) 

C4 = — g[2ai — a2 + 2as — a4 — as — 4(a6 + ay + ag + ag + aio)], 

Cs = — g[ai + a2 + ag + ay — 5(as + a4 + as + ag + ag + aio)]. 

To eliminate the secularity in bracket, we are forced to choose 

2 

Ut^ = K^iu) := -—a^iu), (34) 

so that fl3Tl) finally becomes the following secularity free equation for the first 
correction u^^^: 

M^u^^) = g,. (35) 

This procedure iterates without essential differences at all orders. The terms 
M2U^^^ and (ut^ — K^i^u)) in (^1]) generate, at O^e"^), the terms M2U^^^ and 
M^u^^^ respectively, while the new linear term {ut^ — i{2/6\)uxxxxxx) is re- 
arranged again into the secular factor (w^g — (2/6!)(T6) that must be set to 
zero, to avoid secularities. Since, at 0{e^''), we produce the linear term 
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2 a2k„ 



d'^ u), one infers, in analogy with [T], that u evolves wrt to the 



y^t2k ' •'{2ky.' 

higher times according to the even flows of the NLS hierarchy as follows: 



uu, = K2k{u) :-- 



\k+l 



and one is left with the following triangular set of equations |2], [3] 



M2m(3) = g^, 

M2U^^^ + MiU^^^ = g7, 
Msu^^) + Miu'-^^ + Mqu^^^ = gg, 



(36) 



(37) 



where, f.i., the expression of g^ is presented in formula fl62|) of the Appendix. 
It remains to remark, following [21 E], that the symmetries {c^} and the 
expressions in fl57|) bring naturally to the definition of the following vector 
spaces. 

Definitions Let Vn be the vector space of all differential polynomials in the 
functions {dlu^'^^) and (S^m^*-')) of order n, possessing the gauge symmetry of 
first kind, where: 



order {dlu^^^) = order (9^u('=)) = j + k 



(38) 



(so that order (9^m) = j + 1, since u 



u 



and let Vnijn) be the subspace 



of Vn of all differential polynomials in the functions {diu^''^) and (9^ti(^)), of 
order n, possessing the gauge symmetry of first kind and such that k < m. 
Is is easy to see that, for instance, cT„,ii'„ G P„+i(l), g^ G 7^5(1) and g-/ G 
^7(3) (see il). 



3 Applying the DP integrability test 

Suppose we generate, from the model to be tested, a NLS-type multiscale 
expansion (as in our example); then we have the following scenarios. If such 
a model is S-integrable (C-integrable), 

1) the leading term u of the asymptotic expansion evolves, with respect to 
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the slow times according to the NLS (LS) hierarchy [T]; 

2) there exist elements /n™"* G Vn+m such that the following equations hold 



Therefore equations fj40|) are necessary conditions to be satisfied, in cascade, 
for the model under investigation to be S -(C -) integrable; they are also 
sufficient to guaranty the asymptotic character of the expansion. If equations 
( 140|) are satisfied only up to a certain order, the model under investigation 
is not integrable, being nevertheless "asymptotically integrable up to that 
order" [21 [3]. 

3.1 The DP test and approximate symmetries 

Equations (IH^ . lHOI) . the basic formulae of the DP test, have been derived 
in |2l [3j as a consequence of the existence of a Lax pair for the starting 
integrable model. It follows that, if the conditions fl39|) . p0|) are satisfied up 
to a certain order, the equation under scrutiny admits an approximate Lax 
pair up to that order. 

In this subsection we show how to derive the conditions ([31]), dlOD from 
the existence of infinitely-many symmetries of the starting integrable model. 
This derivation allows one to establish the important relations (to the best 
of our knowledge so far unknown) between the functions /n"^^ G Vm+n of 
the DP test and the symmetries of the starting model. We concentrate our 
attention on the case of difference equations, but our considerations have 
general validity. 

Let ipnt2 ~ ^^ii^n) be an integrable model, say, the AL equation (jl]), and 
let 'ipntirT, ~ J^2m{ipn), m > 2,he one of its infinitely-many higher order com- 
muting fiows (symmetries), reducing, in the continuous limit, to the higher 
commuting flow Ut^^ = K2m{u) of NLS. 

On one hand, from equations (123]) . (1271) . we have that 




+ 



(39) 



(40) 



/ m+k 



) 



'2(m + l) 



u 



{2(m+fe-j)+l) 
'*2j 



(41) 



k>0 \j=m 
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where Ut,^ = K2m{u) (from ([MD) and ^(2|-+*^-i)+i) = + 
j;(2{m+k-j)+i) ^ for some functions f^^^^'^''^^'^^^^ to be specified. On the other 
hand: 

fc>l 

where k!^'^ G P2(m+fc)+i- Equating equations fHTl) and we infer that 
/n"*'' € Vm+ni rn,n & N+ (the basic formula fl39|) of the DP test), and we also 
construct the asymptotic expansion of the generic higher order symmetry 



(m+k-1 / . . , N n\ \ (43) 

j=m ^ ' 

in terms of the NLS higher order symmetries, of their Frechet derivatives 
in the direction of the corrections u^-^-*, j > 1 of the leading term u of the 
expansion fl27j) . and of the output functions G Pm+n of the DP test. 



Therefore, if f^n'^^'^ ^ V2{k+n)+i esists, but fH^^^^ G P2(fc+n)+3 does not, 
Vn G N+, it follows that: 

i) the solution ^(2'^+^) of fl39p is uniformely bounded and the expansion (1271) 
is asymptotic up to the 0{e'^''); therefore the PAE under scrutiny approxi- 
mates well its continuous limit, with an error of O(e^), for time scales up to 
the 0(e-2'=). 

ii) The PAE possesses infinitely-many "approximate" generalized symme- 
tries in the form fH3|) up to the 0(e^^); therefore it is integrable up to that 
order. We remark that, due to the Hamiltonian theory of integrable systems 
[46j-[49j, it is also possible to associate with the PAE infinitely-many "ap- 
proximate" constants of motion in involution, a very useful information in a 
any numerical check. 

3.2 C - and S - integrability at 0{e^) 

In our example, the first of equations fl37|) is already in the form fl39l) . with 
= f!f^ G '^5(1). Assuming now that M^u^'^^ = f['^\ we arrive at the 
consistency 

M4/f = M2/f , (44) 
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that must be viewed as an equation for the unknown . Since = /2 G 
7^5(1), it follows that one must look for /] G Vji^l). The calculation, plain 
but lengthy, has been performed using the algebraic manipulation program 
of Mathematica, and gives the following result. 

Lemma 1. Equation (jUj) admits a unique solution f^^^ G 7^7(1) (presented 
in formula flMl) of the Appendix) iff the coefficients a^-'s appearing in ([2]) 
satisfy the following quadratic constraint 

10 

(ai — 803 — 2a4 — 605 — Sag + Say — Sag — ISag — aio)(^^ aj) = 0. (45) 

Once f[^^ is constructed, f2^^ G ^7(8) (/g^"* = M2W*-^-*) is found from the 
second of equations (1371) : 

if =97- fi'^ (46) 

and is presented in formula fl67|) of the Appendix. 

We first notice the nice factorization of the quadratic constraint fH5l) into 
two linear contraints: 

10 

c = 5]a, = 0, (47) 

i=i 

ai — 3a3 — 2a4 — 6a5 — Sag + Bay — Sag — ISag — aio = 0. (48) 
Therefore we have the following two different scenarios. 

1) If the first constraint p7|) is satisfied by the coefficients a^-'s, the continuous 
limits of dNLS ([2]) is the linear Schrodinger equation. It follows that, in this 
case, equation ([2]) is "asymptotically C-integrable" at O(e^) and one expects 
that, for generic initial data and at time scales of 0(e"^), the dynamics 
according to fl2]). flTr|) be well approximated by the dynamics according to the 
linear Schrodinger equation (l23l) with an error of O(e^). 

In particular, the dNLS ([2]), ([7]) is "asymptotically C-integrable" at O(e^) 

iff: 

ai + 4a3 + 3a6 = 0. (49) 

2) If, instead, the second constraint fHSj) is satisfied by the coefficients aj's, 
the dNLS equation ([2]) is "asymptotically S-integrable" at O(e^) and one ex- 
pects that, for generic initial data and at time scales of 0(e~^), the dynamics 
according to the dNLS equation ([2]), psj) approximates well the dynamics 
according to the NLS equation (OP), ([3]) with an error of O(e^). 
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In particular, i) the dNLS ([2]),® is "asymptotically S-integrable" at O(e^) 
iff the following additional constraint is satisfied: 

as + 2a8 = 0; (50) 

ii) the dNLS (E]),© is "asymptotically S-integrable" at 0{e^) iff the following 
additional constraint is satisfied: 

ai — 8a3 — Qae = 0, (51) 

while the dNLS ([2]), ([8]) is not "asymptotically S-integrable" at 0(6^) (there- 
fore it is not integrable). 

In addition, since the dNLS equation ([2]) is the linear combination of ten 
different discretizations of NLS, it is immediate to check if some of these 
ten discretizations satisfy the constraint fH5]l . Calling dNLS^ the single dis- 
cretization of NLS obtained choosing in ([2]) aj = a^Sjk, j = 1, . . . , 10, it is 
straightforward to see (since the coefficient 02 is the only one absent in (HHIl ) 
that only the dNLS2 equation (coinciding with the AL equation (jlj)) satisfies 
the constraint 0481) (as it has to be, being an integrable system). All the other 
dNLSfc, k ^ 2 equations, including the dNLSi equation do not satisfy 
the constraint fj^Hj) : therefore they are not "asymptotically S-integrable" at 
O(e^) (consequently, they are not integrable) and, for generic initial data 
and at time scales of 0(e^^), their dynamics are expected to be quite differ- 
ent from that of NLS ([I]), ([3]), presumably exhibiting numerical evidence of 
nonintegrability and/or chaos. 

We finally infer that the discretizations ([2]),([H]) and ([2]), ([7]) satisfying re- 
spectively the constraints (!50|l and (ISTI) . the AL equation and any other dNLS 
equation ([2]) satisfying the constraint (148|) are all close to NLS (once the free 
coefficients of each model are normalized to satisfy ([3])) and are all close 
together at time scales of 0(e~^), in the sense mentioned in the introduction. 

It is interesting now to push the integrability test to the next order. Due 
to the above factorization of the constraint fHSl) . the test bifurcates and, in the 
next two subsections, we explore both cases. Before doing that, we observe 
that, given f2^^ G ^5(1) and assuning that the constraint fH5|) be satified, 
equation Mg/i^'* = ^2/6^"* admits a unique solution /g^-* = Mew*-^-' G ^9(1), 
presented in formula (1^91) of the Appendix, and no additional constraint 
appears in this derivation, as predicted by the DP test. 
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3.3 C - integrability at 0{e^) 

Let us assume that the constraint fHTj) be satified. For the construction of 
ff^ = M^u^^^ G ^9(3) from the equation 



we have the following result. 

Lemma 2 Equation fl52|) admits a unique solution ff^ G ^9(3), presented in 
formula fTfTl) of the Appendix, iff the coefficients a/s satisfy the four linear 
constraints ffTTj) . defining a 6-parameter family (but one of these parame- 
ters can always be rescaled away) of dNLS equations ([2]) "asymptotically 
C-integrable" at O(e^). Therefore one expects that, for generic initial data 
and at time scales of 0(e~^), the dynamics according to fl21). f[TT]) well approx- 
imate the dynamics according to the linear Schrodinger equation (!23|) with 
an error of O(e^). 

For instance, the discretization ([21), © satisfies the constraints (fTTj) iff 

Og = — CI3 = (X\. 

The 6-parameter family of dNLS equations (I2]). (ITT]) (or at least some 
particular case of it), being C-integrable at such a high order, is a natural 
candidate to be a C-integrable discrete system. 



Let us assume that the constraint (l48l) be satified. For the construction of 
a unique ff^ = M^u^^^ G "^9(3) from equation fl52l) . we have the following 
result. 

Lemma 3. If the constraint pHj) is satified, equation (1521) admits a unique 
solution /]^^ G 7^9(3), presented in formula flTTj) of the Appendix, iff the 
following five quadratic constraints are satisfied: 



M4/f = M2ft 



(52) 



3.4 S - integrability at 0{e^) 



Qj = 0' J = 1, • • • ,5, 



(53) 



where the Q/s are the following quadratic forms in the 9 variables 
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02, ... , Oiq: 

Qi = — 4a^o + ciioO'2 + 2aioa3 — 0203 + 2a| — 01004 — 202O4 + O3O4+ 
3010O5 — 202O5 — 303O5 — 804O5 — 805 + ISoioOg + 603O6 — 6010O7+ 
402O7 + 603O7 + 404O7 + 2O05O7 + 2406O7 — 8oy + 12oio08 — 302O8+ 
603O8 + 80408 — 905O8 — 606O8 + I807O8 + 2O010O9 — 302O9 — 203O9— 
I304O9 — 2505O9 — 606O9 + 5O07O9 — 24o809 — 24o|, 

(54) 

Q2 = 14oio + 6010O2 + 44010O3 + 402O3 + 26o| + 36010O4 + 502O4+ 
400304 + I704 + 72010O5 + 702O5 + 8803O5 + 6804O5 + 7505+ 
64oio06 + IO02O6 + 7203O6 + 6O04O6 + I2805O6 + 60og — 24010O7— 
202O7 — 2403O7 — I604O7 — 4405O7 — 3206O7 + 4of + 64oio08 + 802O8+ 
6O03O8 + 5404O8 + IO605O8 + IOO06O8 — 2O07O8 + 42o| + I68010O9+ 
2802O9 + 22O03O9 + I7O04O9 + 38205O9 + 33206O9 — IO807O9+ 
2840809 + 466o^, 

(55) 

Q3 = 20o^Q + I5010O2 + 38010O3 — 502O3 — 2203 + 39010O4 — 02O4— 
2303O4 — o| + 63oio05 — II02O5 — 8303O5 — 5204O5 — 75o|+ 
7O010O6 — 14o206 — 7803O6 — 4804O6 — I4805O6 — 84o| — I8010O7+ 
IO02O7 + 4203O7 + 3204O7 + 7605O7 + 8806O7 — 2O07 + 88010O8— 
702O8 — 5403O8 — I504O8 — II905O8 — 1340008 + 8207O8 — 36og+ 
72010O9 — 5902O9 — 30203O9 — 2II04O9 — 53905O9 — 52606O9 + 23407O9 — 
47208O9 - 78809, 

(56) 

Q4 — — 32oiQ — 24oio02 — 56010O3 + 602O3 + 36o| — 7O010O4 — 02O4+ 
3O03O4 + O4 — II4010O5 + O2O5 + 7803O5 + 2O04O5 + 2705 — I2O010O6+ 
I402O6 + 8803O6 + 4804O6 + 8405O6 + 84og + 24010O7 — I402O7 — 6403O7— 
5204O7 — 8O05O7 — II206O7 + 2807 — 164oio08 + 202O8 + 4803O8 + I204O8+ 
I605O8 + 124o608 — II607O8 + 36o| — 22O010O9 + 2202O9 + 2O803O9 + 9604O9+ 
I9605O9 + 29206O9 - 20407O9 + I7608O9 + 300o|, 

(57) 

Q5 = 4o^Q + 3010O2 — 2010O3 — 02O3 — I403 + 3010O4 — 02O4 — I903O4— 

5O4 — 5O10O5 — 3O2O5 — 47O3O5 — 36O4O5 — 39O5 — 20io06 — 6O2O6 — 

3803O6 — 3204O6 — 6805O6 — 36og + 6010O7 + 202O7 + I803O7 + I604O7+ 
2805O7 + 24o607 — 407 + 8010O8 — 302O8 — 3O03O8 — I904O8 — 5905O8— 
62o608 + 2607O8 — 20og — 4O010O9 — 2302O9 — I5O03O9 — II904O9 — 25505O9— 
23O06O9 + 8207O9 — 2I608O9 — 35609. 

(58) 
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The five homogeneous quadratic constraints for nine unknowns, 

characterizing the intersection of 5 quadrics in the real projective space of 
dimension 8, define, in principle, a 4 - parameter family of solutions (but 
one of these parameters can always be rescaled away) whose parametrization 
does not appear to be expressible, in general, in terms of elementary func- 
tions. The corresponding dNLS equation ([2]) is asymptotically S-integrable 
at O(e^) and should well approximate the NLS equation for times up to the 

We observe that, in all these quadratic constraints, 02 is the only coef- 
ficient appearing always multiplied by other coefficients (the term (02)^ is 
absent); therefore the choice 

«j = «2^j2, j = 1, . . . , 10, (59) 

corresponding to the AL equation (jlj), satisfies all constraints, as it has to 
be. Other less trivial explicit solutions of can also be 

constructed, corresponding to the case in which all quadrics degenerate into 
pairs of hyperplanes. Here we display the following two examples: 



«! = — 4a6, a2 = 03 = 4a6, 04 = 0, 
as = -4a6, a? = -^e, ag = ag = aio = 0, 



ai = — 24a9, 02 = 03 = 0, = = —Sag, 

Ctg = lOctg, fly = — Scig, Ctg = — 7cig, CIj^q = 6ctg, 



(60) 



(61) 



corresponding, respectively, to the dNLS equations (fT2l) and (fT3!) presented 
in the Introduction, "asymptotically S-integrable" at O(e^). Therefore one 
expects that, for generic initial data and at time scales of 0(e~'^), the dy- 
namics according to equations (IT^ and (IT^ are good approximations of the 
dynamics according to the NLS equation ([T]), with an error of O(e^), at time 
scales of 0(e~^). Of course, these distinguished equations, passing the test 
at such a high order, are also good candidates to be S-integrable difference 
equations. 

We finally observe that there is no choice of parameters for which the 
dNLS equations ([2]),([H]) and ([2]), ([7]) satisfy the above constraints; therefore 
these two models are not S-integrable at this order (they are not S-integrable 
at all) and do not approximate well NLS at time scales of 0{e~^). 
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4 Summary of the results and future perspec- 
tives 



In this paper we have proposed an algorithmic procedure allowing one: i) to 
study the distance between an integrable PDE and any PAE discretizing it, 
in the small lattice spacing e regime; ii) to test the (asymptotic) integrability 
properties of such a PAE, and iii) to construct infinitely many (approximate) 
symmetries and conserved quantities for it. This method should provide, 
in particular, useful and concrete informations on how good is a numerical 
scheme used to integrate a given integrable PDE. 

The procedure we have proposed, illustrated on the basic prototype ex- 
ample of the nonlinear Schrodinger equation ([1]) and of its discretization ([2]), 
consists in the following three steps: i) the construction of the multiscale 
expansion of a generic solution of the dNLS ([2]) at all orders in e, follow- 
ing pj; ii) the application, to such expansion, of the DP integrability test 
[21 E]; iii) the use the main output of such test to construct infinitely many 
approximate symmetries of the dNLS equation ([2]), through novel formulas 
presented in this paper. 

This approach allows one to study the distance between the integrable 
PDE and any PAE discretizing it. Suppose, for instance, that the asymptotic 
expansion we construct reads ip = u + 0(e"), a > 0, where is a. generic 
solution of the dNLS ([2]) and u is the corresponding solution of ([1]), then if 
the DP test is passed at O(e^), we conclude that: i) the dynamics according 
to the NLS equation ([1]) is well approximated (with an error of 0(e")) by 
the dynamics according to its discretization ([2]), for time scales of 0{t~^); ii) 
the dNLS equation is asymptotically integrable up to that order, construct- 
ing its infinitely many approximate symmetries and constants of motion in 
involution. On the contrary, if the DP test is not passed at that order, the 
dNLS equation is not integrable and one should expect, at the corresponding 
time scale, numerical evidence of nonintegrability and/or chaos. 

We have carried the above procedure up to the O(e^) and we have been 
able to isolate the constraints on the coefficients of the dNLS equation ([2]) 
allowing one to pass the test at that order, in both scenarios of S- and C- 
integrability. 

Numerical experiments to test such theoretical findings are presently un- 
der investigation; preliminary results seem to confirm the theoretical predic- 
tions contained in this paper [JT]. 
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With the same methodology and goals, we are presently investigating 
families of discretizations of the Korteweg-de Vries and Burgers equations 
[52] . other two basic integrable models of natural phenomena. Of course we 
also plan to investigate discretizations of integrable PDEs in which also the 
time variable is discretized. 



5 Appendix 

In this Appendix we display, for completeness, the long outputs of the DP test, obtained 
using the algebraic manipulation program of Mathematica. 
The differential polynomial g-j in ([57]) reads: 

57 = i(liu\u\'^ + l2\u\^u^^^ + ^sUM^^)' + hu^\u\^u^^^ + hu\u^^^\^ + 

Iqu\u\^uI + hulu'-^'^ + lsu\u\'^\Ux\'^ + lg\Ux\'^u'^^^ + lioU^ul + 

heUUxxU^^^ + hllUxl^Uxx + hsUul^ + ligU^\u\'^Uxx + hoUUxxU^^'' + 
hiuluxx + h2u\Uxx\'^ + l23\u\'^uil + hiu'^U^xx + l25UUxUxxx + 
h&UUxUxxx + hlUUxUxxx + h?,\u\'^Uxxxx + h^U^U^x: 



(62) 



where 

h = ~]^c^: h = — fc^i ^3 = 2c, Z4 = I5 = 4c, = 7lio, Ij = ^hi, 

Is = 6Z10, I9 = — + a2 + — 604 + la^ + qq + — llag + 13ag + aio), 
^10 = "M*^^' ^11 = + ^2 + 03 + 04 + 05 - 3a6 + 5a7 - 3(a8 + ag) + 5aio), 
^12 = ^13, ^13 = ^^("1 + a2 + 03 - 5^4 + 7a5 + ag + a7 - llag + 13ag + aio), 

/l4 = 5/10, ^15 = /16, 

/16 = -i(2ai — a2 + 2a3 — a4 - — 4(a6 + a7 + ag + ag + aio), 
li7 = -^(5(ai + a2 + as + a4 + as + ag + a7 + aio) - STag + 77ag), 
^18 = + a2 + as + a4 + as - 8(a6 + a7 + ag + ag + aio)), 

^ig = 2/10, ^20 = -|(ai + a2 + ae + a? - 5(a3 + a4 + as + ag + ag + aio)), (63) 

hi — -jsiO'i + a2 + 03 + 04 + as + ag + a7 - 17(a8 + a9) + 19aio), 

'22 = -Tgo (ll(ai + a2 + as) - 79(a4 + as) + ll(a6 + a7) - 169(ag + ag + aio)), 

^23 = -3 (2ai - a2 + 2a3 - a4 - as - 4(a6 + a7 + ag + ag + aio)), ^24 = 1^20, 

hb = -i^(ai + a2 + a3 + a4 + as + 9a7 - 7(a6 + ag + ag) + 9aio), 

he = "^('^("i + 02 + a3 + aio) — 17a4 + 23as + 3ag + 3a7 — 37ag + 43ag), 

hj = -3^(01 + 02 + a3 — 14a4 + IBas + ag + a7 - 29ag + 31ag + aio), 

^28 = ~ ^"^2 + 2a3 - 3a4 - 305 - 8(ag + a7 + ag + ag + aio)), 

hg = - 180 ("1 + 02 + ag + a7 - 14(a3 + a4 + as + ag + ag + aio)). 

The solution /f ^ e ^7(1) of Af4/f ^ M2.f!t\ where f^^^ = 55 is given in exists 
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unique and reads 
/i^'' = i\aiu\u\ 



aeu^u 



a7Uxxxx\U 
2 



a2Uxx 
2 



OL^Uxxxx^ 
■ OliiUxx\Ux\ 



agUxxxUxU + aioUxxxUxU + aiiUxxxUxU+ (64) 



^Ib^xx^x 



where 



3-(2c2 - C3 + C4 + Scs), a2 = f (4c2 - 2C3 + 6C4 + Scs), 
^(2c2 - C3 + 3c4 + lOcs), a4 = ^(40c2 - IIC3 + 25c4 + 2OC5), 
Y^(8c2 + 5c3 + 7c4 + I6C5), ae = 5|(4c2 + C3 + C4 + Scg), 
^, as = f|, ag = Y2(4c2 + 3c4), aio = y2 (C3 + 2C5), 
an = j^(2c3 + C4), ai2 = j^(3c2 + 2C4), 013=^^(03 + 04 + 405), 

"14 = 12 (2C2 + 5C3 + C4), ai5 = J^(C2 + C3 + 3C5), 



ai 
as 



iff the foUowing constraint is satisfied 



(65) 



2ci - c(2c2 - C3 + c4 + 4c5) = 
on the coefficients Cj's defined in (I33p . This constraint is equivalent to ([451 



(5) 



97 — f!i^^ consequently reads, from (162]) and ([64] 



de-ul-UpM^ + drulu'^^^ + dsulMplu^P + dg\Ux\'^U^^^ + digU^U^ + diiMU2:ui^^ + 



di2UUxUx + disUUxUi; + (ii4|upUj:a; + disUM 

+ d20UUxxu'^^^ + d2lulUxx + C?22w| 
+ C?26'«W x^xxx ~^ d27UUxUxxx'^ 

d2S,\u\^Uxxxx + d2gu'^Uxx: 



(66) 



(67) 
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di = ^(5ai + 2a2 — 4a3 — 04 — ISa^ — ISa^ + liar — Wag — SAag + 2oio), 

d2 = — §c^, ^3 = 2c, di = — c^, ^5 = 4c, 

= ^(95ai + 20a2 + 35a3 + 26a4 - lOGas - 295a6 + 185a7 - 223a8 - 487a9+ 
125aio), dr = -^(ai + 02 + 03 + 04 + 05 - Sae + Sot - 3(a8 + 09) + 5aio), 
ds = ^(llai + 4a2 — 803 — 22a4 — 2a5 — I9ae + 1807 — 55a8 — ISog — 3aio), 
dg = di3, dio = ^(llai + 802 — 13a3 — 22a4 — 10a5 — IQae + 29a7 — 55a8— 
Slag + 5aio), (in = 2^7, di2 = £^13, 

rfi3 = — + a2 + 03 — 5a4 + 7a<^ + ag + 07 — lla8 + ISag + aio), 

di4 = ^(16ai — 20,2 + 13 — 5a4 — 29a5 — 44a6 + 4a7 — 35a8 — 83ag — llaio), 

di5 = die, die = -^(2ai - 02 + 2a3 - 04 - 05 - 4(o6 + 07 + 08 + og + aio)), 

(ii7 = -^{5ai + 2a2 + 5a^ — 2804 + 32a5 — ISag + llaj — ag — 25ag + llaio), 

dis = ^(13ai +a2 + 13a3 + 04 + 05 - llog + 61aj - ll{as + ag) + 61aio), 

dig = ^(lloi + 2o2 — 19a3 — 22a4 — 34a5 — 19ae + ba^ — STog — 61ag — 25aio), 

d20 = —■^(o-i + 02 + flg + 07 — 5(0,3 + O4 + 05 + 08 + 09 + flio)), 

d2i = 55(2ai + 2a2 — 7a3 — 13a4 — 05 — 4a6 + 807 + llag + 35ag — 37aio), 
d22 = i§o(-'-^'^i ~ 0-2 — 4603 — 04 + 59a5 — 1606 ~ 1607 + 44a8 + 164ag+ 
104aio), 

d23 = -^(2ai - 02 + 203 - 04 - 05 - 4(a6 + 07 + 08 + ag + Oio)), ^24 = ^^20, 

d25 = |:(ai + 03 - 06 - 07 - 08 - Og - Olo), 

c?26 = tIo (-'^-'^'^1 ^ ^'^2 + 1103 - 4a4 - 405 - 19(05 + 07 + Og + O9 + Oio)), 

1^27 = Joif^l +0-2 - 403 + O4 - 905 +06 +07 + 6o8 - 14og - 4oio), 
C^28 = TM (I'^'^l - O2 + 1403 - 04 - 05 - 16(06 + 07 + Og + Og + Oio)), 
^29 = lio- 

The unique solution /f ^ = Mqu^^^ e Pg(l) of equation M6/2^^^ = M2/f ^ reads: 

/f ^ = i(/3i|u|8u + (32U,,\uf + /33Wxx|m|^m' + /34M^|U|^M + P5\u^\^\u\^u+ 
Peulu^\u\'^ + ^7Uxa;xa;|u|^ + Ps'U'xxxUx\u\'^U + /3gWa;a:a;Ua: + PloulJu\'^U+ 

Pll\Uxx\'^\u\'^U + /3l2Ua;a;a;U2:|MpW + Pl^UxxxUxU^ + /3l4Wa;a;W^M^ + 
/3l5Ma;a;|Ma;P|Mp + PiqUxxuIu'^ + /3i77;;,' (I,,,?/, + /3i8 jUj; I"*?! + /3l9 (m.t.t + 
/320Ma:a:|Ma:pM^ + /32lMa:xM^ I^P + P22Uxxxxxx\u\'^ + P23Uxxxxxxu'^ + 
P24Uxxxxx\'u\ + P25UxxxxxUUx + (326^ xxxx\'^x\ H~ 027'^xxxx'^xx'^~^ 
I^2SUtxxxUxxU + f329UxxxxUl + P^O^xxxx^xxU + 7/^..,..^.77 + f3z2UxxxUxx^x + 
PsS^xxx^xx^x H" /^34 1'^xa^a: I ^ H" /^Ssl^xa;! '^xx H" PsQ^xxxxx^x^^ 
PsT'^xx'^xxx'^x XXXX 
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where 



/3l = f|(6c2 - 3C3 + 3C4 + 8C5), /?2 = §|(10C2 - 5C3 + 10C4 + I2C5), 

/33 = ig(4C2 - 2C3 + 4C4 + 9C5). /94 = f8(13c2 - 5C3 + 9C4 + IOC5), 
135 = f6(18c2 - 3C3 + 14C4 + 28C5), I3q = f^(5c2 - C3 + 3C4 + Scg), 

/37 = Tfo(6c2 - 3c3 + 15c4 + Tcg), /Sr = 3|o(108c2 - 29c3 + I26C4 + SGcg), 
/Jg = 3fo(36c2 + 7c3 + 66c4 + 52C5), /3io = 7|j(152c2 - 41c3 + 169c4 + 84c5), 
/3n = 3|o(44c2 - 7c3 + 79c4 + II8C5), /3i2 = 3fo(16c2 + 7c3 + 26c4 + 62C5), 

/^13 = 3I0 (8C2 + C3 + 6C4 + 26C5), /3l4 = 75(33C2 - 7C3 + 21C4 + 14C5), 

/3i5 = 3fo(278c2 + llc3 + 276c4 + 236C5), /3i6 = 3fo(39c2 + 8C3 + 33c4 + 68C5), 

/3l7 = ^(22C2 - C3 + 13C4 + I6C5), /3l8 = 7§o(158C2 + 31C3 + IOIC4 + 176C5), 

/3l9 = tIo (12C2 - C3 + 9C4 + 44C5), /320 = lio (16C2 + 7C3 + 12C4 + 42C5), 

/321 = 3|J(114C2-12C3 + 103C4 + 158C5), ^22^^, 1^23 = jk^ 

/324 = 2io(2c3 + 3c4), /325 = 3Ho(f^3 +4C5), /?26 = tM (18C2 + 21C3 + 25C4) , 

/327 = 3go(15^2 + I3C4), /?28 = 75o(9c3 + llC4 + I2C5), 

/?29 = 36o(C2 + 2C3 + IOC5), = 3^(2C3 + C4 + IIC5), 

/331 = Ti4(4C2 + 3C4), /332 = 755 (50c2 + 35c3 + 34C4), 

/?33 = 3|o(llC2 + 17C3 + 10C4 + IScg), /334 = 750(11^3 + 4C4 + I8C5), 

035 = 7io(20C2 + 25C3 + IIC4 + 2OC5), 

036 = 2i0 (4C2 + 5C4), /337 = 75o(8C2 + 31C3 + 4C4 + SOCs), 
/^38 = sio (2C2 - C3 + 5C4 + I4C5), 

and no additional constraint on the coefficients Cj's appears. 
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The unique solution /i^^ = M^u^^^ G ^9(3) of equation Mi/f ^ = M2/4^^ reads 

/i^^ = i(Siu\u\^ + d2Uxx\u\'^ + 63Uxxu'^\u\'^ + diululul"^ + 65\Ux:\^u\u\'^+ 
6eul\u\'^u^ + S7Uxxxx\u\'^ + SsUxxxUx\u\'^u + 6qUxxxUx\u\'^u + dioulJu\'^u+ 

Sll\Uxx\'^\u\'^U + di2UxxxUx\u\'^U + Si^UxxxUxU^ + 5uUxxUlv? + 5i5Uxx\Ux\'^\u\'^ + 
heUxxUiu^ + SnulUxU + Sis\Ux\'^U + Sigul^U^ + S20Uxx\Ux\'^u'^ + S2lUxxUl\u\'^ + 
^22'^xxxxxx\'^\ ~^ ^23'^xxxxxx'^ "1" "I" ^25^xxxxx^x'^ ~t~ ^26^xxxx I'^x | "I" 

^27^xxxx^xx'^ ^28^xxxx'^xx^ xxxx 

^32'^xxx^xx^x ~l~ ^33^xxx^xx^x H" ^34|Wa;a;a;| ^ H" *^35^a:a: |^a:a: | H" ^36^xxxxx'^x^~^ 

hlUxxUxxxUx + (^SSM^IwpWxxxx + 7l^txxxx|up + l2UxxxUxU + 73MxxxMxW+ 

(3) - , (3) I 12 , (3) - , (3) - , (3) _ , 

liUxxUxxU + 75WXX \Ux\ + ^eUxxUUxx + 77Wx ^xxx" + 78'Wx 'WxxWx + 

^gU^x^UxUxx + IwU^x^UUxxx + 7llW^^^UxxxxW + 712^^^^ "xxx^x + 713"'^^^ |Uxx|^ + 

^lAtl'^^^llxUxxx + 715w''^-'mMxxxx + 7l6uixxxU^ + 7l7uixxMxU + 7l8ulx MxxU+ 

7l9UxxUx + 720Ux^''mxxxU + 72lUx^^MxxUx + l22U^^^ UxxxxU + J23U^^^ UxxxUx + 

724M^^^mL + 725WXX l^l* + 726«x^^Wx|m|^M + 727Mx^^Wx 

728u(^^Uxx|upU + 729U^^^Mx'"^ + ISOU^^Auxl"^ H'^ + 73lW^^^UxW^ + 732M^^^Mxx 
733wix IwP'^^ + 734Mi^''Mx|upU + JSbU^x^ UxU^ + JsqU^^'^ Uxx\u\'^U + 737^^^) + 
738U^^VxpM^ + 739U^^^MxxU^ + 740U^^Vl® + l4lU^^'> \u\^U^ + aiU^xxU^^'>U+ 
a2ui^ U + azui\^'^'^Ux + Cr4w(^) Uxx + CTsUxx W^^^U + (Teluif^ |^w+ 
C77uif'«(^)«a; + crg'U^^^'uix^U + CTgU^^^uif '«x + crio I^U^^a; + CTn 
C7i2|u(=^)|2|u|2u + <713(«(3))V), 

(71) 
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where: 

51 = gfg (336rfi + 4(6c2 - 3c3 + 2c4 + 1005)^2 + 605(202 - C3 + 2c4 + 405)^3+ 

4(6c2 - 3c3 + 3c4 + 2c5)d4 + (-4c| + 2C2C5 - C3C5 + C4C5 - 2c^)d5+ 
c(2(14c2 - 5c3 + 9c4 + 28c5)dii - 2(10c2 - 3c3 + 7c4 + 8c5)c;i2 - 4(6C2- 

3C3 + C4 + 16C5)rfi3 + 16di4 + 2(18C2 - 7C3 + llC4 + 30C5)rfi5 - 2(6C2- 

3c3 + 3c4 + 14c5)rfi6 - 48di9 - 4(6c2 - 3c3 + 3c4 + 805)^20 + 4(-4c2 + C3+ 
2c4 + 15c5)d23 + 4(-6c2 + 5c3 + C4 - 605)^24 + 32^ - 164 + 4(2c2 - C3+ 

C4 + 4:C5)dg) + 8c2(-2(ii7 + 4(ii8 + 2^22 + 7^25 - 7^26 - 5^28 + 16^29)) , 

52 = (96di + 10c4(i2 + 4c4C5(i3 - 6c5(i4 - 2c|(i5 + c(2(4c2 - C3 + 6C4+ 

8c5)dii - 2(2c2 + 5ci)di2 - 2(6c2 - 3c3 + C4 + 14c5)(ii3 + 104di4 + (18c2- 

7c3 + IIC4 + 24c5)di5 + 3(-2c2 + cg - C4 - 4c5)di6 - 24di9 + (34c2 - 19c3+ 
31c4 + 80c5)d23 + 2(-6c2 + 5c3 - 5c4 - 10c5)rf24 + 164 - 8ds + 2(2c2- 
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C3 + C4 + 405)4) + 4c^(-2di7 + 4(ii8 + 242 + 1345 - 1346 + 48 + 4d: 
^3 = ^(12di + 2c54 + 2c44 + c(2c5dii - 2c542 + Sdii + 24di9 + (2c2- 
C3 + C4 + 4c5)4o + 2(2c2 - C3 + C4 + 5c5)43 + 5(2c2 - C3 + 04)44)+ 
8c2(45 - 46 + 48 - 249)), 

'^4=124 (215^1 + 2(3C2 + 6C4 - 2C5)4 + 2(2C3C4 + 2C2C5 + C3C5 + C4C5- 

2ci)4 - 24c54 - (6cl - 3C3C5 + C4C5 + 2ci)4 + c(-244o + (52c2 - 19c3+ 

9c4 + 54c5)dii + 2(-7c2 + 3c3 - 8C4 + 9c5)42 + 2(-10c2 + 4c3 - 5c4- 

22c5)43 + 2844 + 2(33c2 - 4c3 + 5c4 + 16c5)cii5 + 2(-22c2 + 12c3 - I8C4- 

19c5)di6 - 102di9 + (4c2 - IOC3 - 6C4 + 7c5)4o + (154c2 - 58c3 + 121c4+ 

199c5)43 + 3(-22c2 + I8C3 - 3c4 - 22c5)44 + 1324 + (lOca - 5c3 + 5c4+ 

4c5)4 - 164 + 2(6c2 - 3c3 + 3c4 + 13c5)4) + 2c^{-A2dn + 118(ii8 - 84i+ (72) 

1742 + 5445 - 3746 + 747 + 19348 + I6O49)) , 

^5 = ^ (48di + 2(3c3 - 2c4 + 4c5)4 + 2C4C54 + 2(-C3 - 6C4 + 4c5)4+ 

2(c| - ci)4 + c(-84o + 4(c2 + 4c5)4i + (IOC2 - 7c3 + IIC4 + 805)42+ 
(6c2 - 3c3 + 3c4 + 4c5)di3 + 24rfi4 - 4(3c2 - C3 + C4 + 05)^15 + 4(4c2- 
2c3 + 3c4 + 3c5)dl6 + 16di9 + 4c54o + 2(28c2 - 15c3 + 6C4 + 42c5)43+ 

4(-4c2 + C3 + 5C4 - 17C5)44 + 84 + 324 + (-2C2 + C3 - C4 - 4C5)4) + 

4c2(2di7 + 4rfi8 - 642 + 345 + 46 - 447 - 1348 - 449)) , 

S6 = jg (124 + 2ci4 + 2(c2 - C4)4 - C4C5d5 + c(204o + 2(2c2- 

C3 + 4C5)42 - 4rfi4 + 2(C2 + 2c5)di5 - 2(C4 + C5)46 + 8di9 + (8C2- 
2C3 + 7C4 + 8C5)43 + 2(10C2 - 5C3 + llC4 + 5C5)44 + 44 + 44) + 

2c2(-2di7 + 12di8 - 242 + 45 + 346 - 447 + 948 + 1249)) , 

(57 = ^ + 2C443 - C544 + 2c(45 - 46 + 948 - 249)) , 

^8 = ^ (cidii - C5rfi3 + 844 + 2C245 + (6C2 + 7C4)43 - 2C544+ 

84 + 2c(-2di7 + 8di8 + 1245 - 246 - 47 + 2748 - 449)) , 

^3 = l4 ("^4^12 + 4(ii4 + 2c5di5 - 2C446 + (3C3 + 4C4)43 " C344 + 44+ 

2c(448 - 242 + 45 + 946 - 47 + 1848)) , 

^10 = 53 (204^11 - 2c5di3 + 8rfi4 + 3c4ii^- c^die + 2c243 + C443 + 64+ 
2c(ll48 - 42 + 245 - 247 + 648)) , 



^11 24 

11 - 2c4rfi3 + 4di4 + 05^15 + 04^16 + 12rfi9 + 04^20 + (4C3 + C4+ 
6C5)d23 + (-C4 + 2C5)d24 + 2^ + H^IS + 4^22 + ^25 - ^26 + 13^28 - 2(^29)) , 
S12 = 53 (^Csrfii + C4(ii3 + 6rfi9 + 03^23 + 1005^23 + £3(^24 - 404^24 + 2^8 + 
2c(d25 + £^26 + 7^27 + 8(^28 + 8(^29)) , 
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5j ^4^10 + 05^12 + 2c?i9 + 202^24 + 204^24 + 2c(d25 + £^26 + £^27 + 8(^29)) , 

Si4 = ^ (4(3c2 + C4)dii + 8rfi4 + 2(2c2 - 04)^15 - 2c5die + (12c2 + 704)^23+ 

2c5(i24 + 244 - 605^7 + 2c(-4di7 + 20di8 - 6^21 - 2^22 + 37^25 - 2^26 + 

37^28 + 4^29)) , 

'^15 = ii ((5C3 + f-4)dll + (2C2 + C4)di2 + (-3c3 + 2cr,)dui + 8di4 + (C3 + 
4C5)(ii5 + (C3 - 4c4)cil6 + (8C2 + IIC3 + 6C4 - 4c5)ci23 + 4(-C3 + C4- 
2C5)d24 + 8^6 - 404^7 + 20^ + Cida + 4c(4di7 + 8rfi8 + 11(^25 + 5(^26 - 7^27 + 
34rf28-8rf29)), 

^16 = ^ (l2dio - 2C5rfii + 2(2C3 + C4 + 3C5)rfl2 - 2(2C2 + C4)di3 + 

4di4 + 2c2di6 + 6(ii9 - 305^20 - (2c3 + 705)^23 + (4c2 - 2c3 + C4+ 

4c5)d24 + 8^8 + 2c5d9 + 2c(2di7 + 6di8 + 2d2l + d22 - 3(i25 + 26^26 - 11^27- 
14(^28 + 4^29)), 

Sl7 = ^ ((C2 + 3C3)dii - 2c3di3 + 2(C2 + C3 - 05)^23 + 204^24 + 4(^6 + 
(C3 - 2c4)rf7 + 6^8 + 02^9 + 4c((ii7 + 2(^21 + 7^25 " 4^27 + 6^28)) , 

^18 = ^ (6diQ - 2c5(iii + (c2 + 3c3 + 605)^12 - (3c2 - C3 + 204)^13+ 

6rfl9 - 3C5rf20 - (2C3 + 7c5)rf23 + (2C2 - 6C3 + 04)^24 + 2de + 02^7+ 

8^8 + Cadg + 2c5d9 + 2c(4di7 + 6rfl8 + 2^21 - 3^22 - 3^25 + 20^26 - 5^27 

-14(i28 - 20^29) 




2^ [ 2dio + 2rfi9 + C5(i20 + (2C2 + C4)d24 + 2c(di8 + ^22 + 6(i29)) , 



; 24dio + (c3 + 2c5)(ii2 + 2c2(ii3 + 04^13 + 12rfi9 + (C3+ 

3c5(i23 - (14C2 - 3C3 + lCi)d2i + Ms + Csdg + 2c(2(ii7- 
6rfi8 + 2^21 + 7rf22 - 3rf25 + 6^26 + 9^27 - 36(^29)) , 

^21 = ^ ((C3 + 6C5)dii + (C3 - 3C4)rfi3 + 18di9 + £2^20 + (C2 + 6C3 + 10C5)d23- 
4(C3 + C4)rf24 + 4:de + CAd7 + Ads + 4c(4rf21 + ^22 + 5^25 - 2^26 + 3rf27 + 20^28- 

16d29)), 

<522-^, '^23 = ^, (^24 = 1^(2^26+^28), '525 = 1^(^27 + 2^29), 
^26 = ^(2dl7 + d25 + 3(i26 + C?28), 1527 = ^(4(ii8 + 3^25 + 2^28), 
^28 = 13(2^22 + d26 + d2s), S29 = 1^(^21 + d27 + M29), 
ho = 1^(^22 + ^27 + 4^29), (531 = l^(3dl8 + 2(^25), 

^32 = ]^(5(il7 + 2C?18 + d25 + 2(i26), <533 = -^{dn + Ad2l + 3^22 + C^25 + d2Q), 
hi = ^{d22 + ^26 + 2c?27), <535 = l^(c^l7 + ^18 + 3(^21 + 2(^22), 
'536 = 1^(2(^25 + 3(f28), hi = ]^(C^17 + 2(^21 + ^22 + 5(^27), 
^38 = ^(2(il9 + ':5C'23 + C;4'i24 + 4c((i28 + 3(^29)), 
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71 = ^, 72 = i^(2dii+3d23), 73 = i^(2di2+d23), 

74 = ]^(3(iii + 2di5 + 2(^23), 75 = + 3di2 + ^23 + 2^9), 

76 = ]^((ii2 + 2^20 + £^23), 77 = ^(2*^11 + 3^15), 

78 = ^(liii + 2di2 + di5 + 3dg), 79 = ^(dn + rfi2 + 8^20 + ^9), 

710 = ^(dl2 + d2o), 711 = 712 = ]^(C^15 + 2^9), 

713 = ^{dl5 + 2^20 + dg), 714 = ]^(rf20 + ^g), 715 = if, 

716 = if, 717 = ]^(rfl3 + 2(i24), 718 = 1^(C^13 + cfl6 + 4^24), 

7i9 = T2('^i3 + 3^24 + £^7), 720 = ]^(2di3 + cJie), 721 = 12(5^13 +c^i6 + 2^7), 
722 = 723 = ]^(3(ii6 + 4^7), 724 = ^{2die + idr), 

725 = 1^(2^2 + C(rfii - rfi2 + 9d23 - 2(^24)), 

726 = J2 (^^2 + 204^3 - Csrfs + c(llrfii - di2 - ^13 + 2di5 - 2^20 + 10^23- 
8^24)), 

727 = ]^(2(i2 + 2c5(i3 - C4(i5 + c(dii + 7di2 + ^13 + 2du, - 2diQ + 8^23)), 

728 = 5j (8^2 + 604^3 - 05^5 + 2c(2dii - 2rfi3 + lldi5 - die - d20 + 2^23)) , 

729 = ^ (6d2 + 2(c2 + C3 - 2c5)d3 - (c3 - €4)^5 + 2c(4dii - 2di3 + 4di5- 



2^20 + 2c?23 + d7 + d9)^ 

730 = M (8^2 + 2(c3 + 4c5)d3 + (C3 - 4c4)d5 + 4c(dii + (ii2 + di3 + 4di5- (74) 

2di6 + 2^20 + 6d23 - 4d24 + 4^9)) , 

731 = ^(2^2 + C2d5 + 2c(2di2 + 2d2o + 2d24 + dr + dg)), 

732 = ^(4c/2 + 2c5d^ + Cid^ + 2c(di5 + ^16 + 7^20 + 2(^23 + 2(^24)), 

733 = ^(rf4 + 2c(rf23 + 3(^24)), 

734 = 1^(3^4 + c(dii + di2 + 7rfi3 + 4(^23 + 4(^24)), 

735 = ^{di + C(dii + di2 + Cil3 + 4(^24)), 

736 = ^ ^12^4 + Cidr, + 2c(2rfii - 2rfi3 + Srfis + 7di6 - ^20 + 6^23 - 2^24)) , 

737 = ^ f 18^4 + 02^5 + 4c(3dii - 2rfi2 + rfi3 + rfi6 + 7(^23 - 8^24 + 4^7) , 

738 = s (l8d4 + C3d5 + 4c(d7 + rfg - Bdn + 4di2 + di3 + 2rfi6 - 3^23 - 10^24)) , 

739 = ^{2di + C5d5 + 2c{di5 + die + ^20 + 2^24)), 

740 = ^ (80d2 + 2(6c2 - 3c3 + C4 + 12c5)d3 + 16^4 + (10c2 - 5c3 + 5c4+ 
805)^5 + 8c(3dii - di2 - 2di3 + 3di5 - c^ie - 2d2a + 2^23 - 2^24)) , 

741 = ^ (8^2 + 24(i4 + (2C2 - C3 + C4 + 4c5)(i5 + 8c(dii - di2 + d23 - 2^24)) , 

„ ds „ da „ da „ da „ „ dg „ „ ds 

CTl = -f , 0-2 = ^, 0-3 = -^, 0-4 = jf, cr5 = crio = -f, CTe = (77 = ^, 

(T8=<79 = f|, (711 = ^(7(^3 + ^5), Cri2 = §(^3 + 4c?5), Cri3 = ^((^3 + cig), 

iff the following three sets of constraints are satisfied: 

di5 — diQ = ds — 2d3 = adz — 2cc?23 = — <^i2 + di^ — d\5 + d2o + ^23 — 2c?24 = 

-dg + dii + di3 - 2d7 = (2c2 - C3 + C4)ci3 - 2{dii - du + die - (^20 + 2rf24)c = (75) 

(C2 - C3 + C4 + C5)d3 + 2{di3 - die - (^24 - d7)c = 0, 

(2C2 - C3 + C4 + 2c5)d3 - 2(^4 + M24C = -M2 + (6C2 - 3C3 + 3C4 + 8C5)(i3 + 8c!24C = 0, (76) 
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-4dio + 2c5(ii3 - 2c5rfi6 + 2dig + 05^20 + 05^23 + (2c2 -04- 405)^24+ 

2c(2di8 - d22 ~ d25 + 2^26 - ^27 - 2^28 + 4^29) = 0, 

-4di4 - 2c4(iii + 2(c2 - C4 + 05)^16 + (4c2 + 3ci)d23 + 205^24 + 2c(-2di7+ 

4di8 + 2^22 + ^25 - d26 + 13^28 + 4^29) = 0, 

-2^8 - (c4 - 2c5)dj - 2(c4 - C5)(iii + (2c2 + cg + 2c4 + 05)^13 - (C3 + 804)^16 + 

6^19 + (-C3 + 2C4 + 05)^20 + (5C2 - C3 + 3c4)(i23 - (4C2 + C3 + 6c4)d24 + 
2c(2rfl7 - 2di8 - 4^21 - 7^22 - d25 + 2^26 + 11^27 + 10^28 " 20^29) = 0, 

(77) 

-4d6 + (2C2 + 2C3 -C4 -4c5)«i7 - 12dio + (C3 + 2c5)dll - (C2 -4c4 + 9c5)di3- ^ ' 

4di4 + (3c2 - 5c5)di6 + 12di9 - (2c2 + C3 + 2c4 - 405)^20 + (c2 + 4c3 + C4- 

C5)d23 + (6C2 - 2C3 + 3C4 - 12c5)rf24 + {Sdn + 20dl8 + 4^21 - 14^22 + 10rf25 + 
4^26 + 10^27 + 44^28 - 16d29) = 0, 

-24di + 6{cl - C4Cr,)d3 + 2(6c2 - £4)^4 + (2(2c2 + 2c3 - C4 + 405)^7 - 4^- 
52dio + 2(-4c2 + C3 - C4 + 2c5)dii + 2(c2 + C3 - 6C4 + 605)^13 - 16rfi4 + 2(5c2— 
C3 + C4 - C5)di6 + 62di9 + (-4c2 + 2c3 - 12c4 + 505)^20 + (2c2 + 12c3 + 4c4+ 

13c5)d23 + (6C2 - 12C3 + 5c4)d24)c + 2c2(-16di7 + 50di8 + 4^21 - 13^22 + ^25 + 
20^26 + 3^27 + 52^28 + 4^29) = 0. 

The first set of constraints ([75)1 is automatically satisfied by the parametrizations ([55)) 
and (|68p . while the second set of two constraints ([76)) is satisfied by the parametrizations 
([33)) and ([68| and by the O(e^) constraint ([45)) . The remaining five constraints ([77)) are 
equivalent to the five quadratic constraints ([53| .([54 )) -([58 | in the S-integrability scenario 
in which ([48| holds, and to the linear constraints ([TT)) in the C-integrability scenario in 
which c = (the last constraint is automatically satisfied by the condition c = and, 
in the remaining constraints, the quadrics degenerate into the hyperplanes described by 
equations ([TT)) ). 

Acknowledgements. We acknowledge interesting discussions with U. Agli- 
etti. 



References 

[1] A. Degasperis, S. V. Manakov and P. M. Santini, "Multi-scale pertur- 
bation beyond the nonlinear Schrodinger equation. 1", Physica D 100, 
187-211 (1997). 

[2] A. Degasperis and M. Procesi, 'Asymptotic integrability" in Simmetry 
and perturbation theory, SPT98 23-37, edited by A. Degasperis and G. 
Gaeta, World Scientific, Singapore (1999). 



30 



[3] A. Degasperis, "Multiscale expansion and integrability of dispersive wave 
equations" , lectures given at the Euro Summer School "What is integra- 
bility?", Isaac Newton Institute, Cambridge, U.K., 13-24 August (2001); 
in Integrability edited by A. Mikhailov, Lecture Notes in Physics 767, 
Springer, Berlin- Heidelberg (2009). 

[4] P. Winternitz, "Symmetries of discrete systems. Discrete integrable sys- 
tems", 185-243, Lecture Notes in Phys. 644i eds. B. Grammaticos, Y. 
Kossmann-Schwarzbach and T. Tamizhmani, Springer, Berlin, 2004. 

[5] F. Valiquette and P. Winternitz, "Discretization of partial 
differential equations preserving their physical symmetries", 
arXiv:math-ph/0 50706ll 

[6] T. R. Taha and M. J. Ablowitz, "Analytical and numerical aspects of cer- 
tain nonlinear evolution equations II. Numerical, nonlinear Schrodinger 
equation", J. Comput. Phys. 55, 192 (1984). 

[7] B. M. Herbst and M. J. Ablowitz, "Numerically induced chaos in the 
nonlinear Schrodinger equation", Phys. Rev. Lett. 62, 2065-2068 (1989). 

[8] P. L. Kelley, Phys. Rev. Lett. (1965) 15 1005. 

[9] V. E. Zakharov, Soviet Phys. JETP (1968) 994-998. 

[10] D. J. Benney and A. C. Newell, J. Math, and Phys. (now Stud. Appl. 
Math.) 46 (1967) 133-139. 

[11] A. Hasegawa and T. Tappert, Appl. Phys. Lett. 23 (1972) 142. 

[12] H. Hasimoto and H. Ono, J. Phys. Soc. Japan 33 (1972) 805. 

[13] T. Taniuti, Suppl. Proc. Th. Phys. 55, 1 (1974). 

[14] Y. Kodama and T. Taniuti, J. Phys. Soc. Japan, 45, 298 (1978). 

[15] F. Calogero and W. Echkhaus, "Nonlinear evolution equations, rescal- 
ings, model PDEs and their integrability: I" Inverse Problems 3, 229-262 
(1987). 

[16] V. E. Zakharov and A. S. Shabat, Soviet Phys. JETP 34 (1972) 62. 



31 



[17] V.E. Zakharov and S. V. Manakov, "Resonance interaction of wave pack- 
ets", Soviet Physics JETP 42 (1975). 

[18] V. E. Zakharov and E. A. Kuznetsov, Physica D 18, 455 (1986). 

[19] D. J. Korteweg and F. de Vries, " On the Change of Form of Long Waves 
Advancing in a Rectangular Canal, and on a New Type of Long Station- 
ary Waves." Philos. Mag. 39, 422-443, 1895. 

[20] C. S. Gardner, C. S. Greene, M. D. Kruskal, and R. M. Miura, "Method 
for Solving the Korteweg-de Vries Equation." Phys. Rev. Lett. 19, 1095- 
1097, 1967. 

[21] F. Calogero, "Why are certain nonlinear PDEs both widely applicable 
and integrable?" in What is integrability? 1-62, edited by V.E. Zakharov, 
Springer, Berlin- Heidelberg (1991). 

[22] E. Hopf, Commun. Pure Appl. Math. 3, 201 (1950). J. D. Cole, Quan. 
Appl. Math. 9, 225 (1951). 

[23] V.E. Zakharov, S.V. Manakov, S.P.Novikov and L.P.Pitacvsky, Theory of 
Solitons. The Inverse Problem Method^ Plenum Press (1984) 

[24] F. Calogero and A. Degasperis, "Spectral Transform and solitons: tools 
to solve and investigate nonlinear evolution equations. Volume one". 
North Holland, Amsterdam (1982). 

[25] M.J.Ablowitz and P.C.Clarkson, Solitons, Nonlinear Evolution Equa- 
tions and Inverse Scattering, Cambridge University Press, Cambridge, 
1991. 

[26] B. Konopelchenko, Solitons in Multidimensions, World Scientific, Sin- 
gapore (1993). 

[27] F. Calogero and W. Echkhaus, "Nonlinear evolution equations, rcscal- 
ings, model PDEs and their integrability: II", Inverse Problems 4, 11-13 
(1988). 

[28] F. Calogero, A. Degasperis and X.D. Ji, "Nonlinear Schrodinger-type 
equations from multiscalc reduction of PDEs. I. Systematic derivation" , 
J. Math. Phys. 41, 6399-6443 (2000). 



32 



[29] F. Calogero, A. Degasperis and X.D. Ji, "Nonlinear Schrodinger-type 
equations from multiscale reduction of PDEs. II. Necessary conditions of 
integrability for real PDEs", J. Math. Phys. 42, 2635-2652 (2001). 

[30] Y. Kodama and A. V. Mikhailov, "Obstacles to asymptotic integra- 
bility", in Algebraic aspects of integrable systems: in memory of Irene 
Dorfman, 173-204, edited by A. S. Fokas and I. M. Gel'fand, Birkhauser, 
Boston (1996). 

[31] D. Levi, M. Petrera and C. Scimiterna, "On the integrability of the 
discrete nonlinear Schrodinger equation", E. P. L. 84, 10003 (2008). 

[32] C. Scimiterna, "Multiscale techniques for nonlinear difference equa- 
tions", PhD Thesys, Dept. of Physics, University of Roma3, Roma, Italy 
(2009). 

[33] C. Scimiterna, "Multiscale reduction of discrete Korteweg-de Vries equa- 
tions", J. Phys. A: Math. Theor. (special issue for the Conference SIDE 
8). 

[34] R. I. Yamilov, "Symmetries and integrability criteria for differential dif- 
ference equations". J. Phys. A: Math. Gen., R541-R623 (2006). 

[35] V. V. Sokolov and A. B. Shabat. "Classification of integrable evolution 
equations", Sov. Sci. Rev. section C4, 221-80 (1984). 

[36] A. V. Mikhailov, A. B. Shabat and R. I. Yamilov, "The symmetry ap- 
proach to the classification of nonlinear equations. Complete lists of in- 
tegrable systems", Russian Math. Surveys 42/4, 1-63 (1987). 

[37] D. Pehnovski, "Translationally invariant nonlinear Schrodinger lattices" , 
Nonhnearity 19, 2695-2716 (2006). 

[38] M. J. Ablowitz and J. F. Ladik, J. Math. Phys. 17, 1011 (1979). 

[39] A. S. Davydov, J. Thor. Biol. 38 559 (1973). 

[40] W. P. Su, J. R. Schieffer and A. J. Heeger, Phys. Rev. Lett. 42 698 
(1979). 

[41] J. C. Eilbeck, P. S. Lomdhal, A. C. Scott, Physica D 16, 318 (1985). 



33 



[42] D. Hennig and G. Tsironis, Physics Reports 307, 333 (1999). 

[43] F. K. AbduUaev, B. B. Baizakov, S. A. Darmanyan, V. V. Konotop and 
M. Salerno, Phys. Rev. A 64, 043606 (2001). 

[44] M. Oster, M. Johansson and A. Eriksson, "Enhanced mobihty of 
strongly localized modes in waveguide arrays by inversion of stability", 
Phys. Rev. E 67, 056606 (2003). 

[45] C. Claude, Y. S. Kishar, O. Kluth and K. H. Spatschek, "Moving local- 
ized modes in nonhnear lattices", Phys. Rev. B 47, 14228-14232 (1993). 

[46] F. Magri, J. Math. Phys. 19, 1156 (1978). 

[47] I. Gel'fand and I. Dorfman, Funct. Anal. Appl. 13 (1979); 14 (1980). 

[48] A. S. Fokas and B. Fuchssteiner, Lett. Nuovo Cimento 28 299 (1980); 
Physica 4D 47 (1981). 

[49] P. M. Santini and A. S. Fokas: "Recursion operators and bi-hamiltonian 
structures in multidimensions.I" ; Comm. Math. Phys. 115, 375-419 (1988). 

[50] M. J. Ablowitz, D. Kaup, A. C. Newell and H. Segur, Stud. Appl. Math. 
53, 249 (1974). 

[51] U. Aglietti and P. M. Santini, "Multiscale expansions of difference equa- 
tions in the small lattice spacing regime. Integrability test and numerical 
confirmations", preprint (in preparation). 

[52] U. Aglietti, P. M. Santini and C. Scimiterna, preprint (in preparation). 



34 



